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Abstract. The gaseous giant planets WASP-4b and WASP-5b are transiting 12-magnitude solar-type stars in the 
Southern hemisphere. The aim of the present work is to refine the parameters of these systems using high cadence 
VLT/FORS2 z-band transit photometry and high-resolution VLT/UVES spectroscopy. For WASP-4, the new 
Oh' estimates for the planet radius and mass from a combined analysis of our VLT data with previously published 

2 ■ Pv = °- 55 iao2 PJ- The radius and mass for the host star are R* = 0-$>7to%t Re an d M* = 0.85^°;^ M G . Our 

ground-based photometry reaches 550 ppm at time sampling of ~ 50 seconds. Nevertheless, we also report the 
presence of an instrumental effect on the VLT that degraded our photometry for the WASP-5 observations. This 
effect could be a major problem for similar programs. Our new estimates for the parameters of the WASP-5 
system are R p = 1.09 ±0.07 Rj, M v = 1.58±g ig Mj, p p = 1.23tg;?e PJ, R* = 1.03±g:gf Rq, and M» = 0.96±g;Jl 
Mq. The measured size of WASP-5b agrees well with the basic models of irradiated planets, while WASP-4b is 
qq ■ clearly an 'anomalously' large planet. 

On ' 

, Key words, binaries: eclipsing - planetary systems - stars: individual: WASP-4 - stars: individual: WASP-5 

techniques: photometric - techniques: spectroscopic 

OO ! 1- Introduction Station of the South African Astronomical Observatory, 

■ SuperWASP-South brings a second field of view of 482 
• • So far, the planets that transit their parent stars have gquare degrees tQ the suryeyj allowing it to search for tran _ 

. . ^ undoubtedly brought the most important pieces of mfor- siting planetg in ft large portion of the gky 

■ mation about the physics and composition of the plan- _,, , /,„., , „„„„ , 

. , , . , , f\ y o i o ± / - u The planets WASP-4b Wilson et al. 2008, hereafter 

<-h etary objects outside our Solar System (see review by „ T „„, , , , „„„„ , „ 

CO ™ ■ , , „ n(1 M „ , W08) and WASP-5b Anderson et al. 2008, hereafter A08 

Charbonneau et al. 2007). Most of the transit detections ' ... . . , _ ' 

, , r j-jcu were the nrst transiting planets detected by Super WASP- 

are due to a tew ground- based wide- field surveys target- , ml , , . , 

ing stars brighter than V~13: HAT (Bakos et al. 2004), Soutk Th ^ are both S as S iants sll S htly *f avier * ha * 

SuperWASP (Pollaco et al. 2006), TrES (O'Donovan et Jupi o ter and <f bltm S very cl °f ( P = ^Sand J-628 days) 

i nnnc\ j vr\ /h ri n u j. 1 onnc^ a ,1 to 12 mag solar-type stars. 1 he analysis of the WASP and 

al. 2006), and XO (McCullough et al. 2005 . Among these „ „ b . . , ,. , „ , , , , „ , 

a «ricTi • n 1 • A, i . i follow-up data led to radius values of about 1.4 and 1.1 

surveys, SuperWASP is the one showing the largest har- „ TTT . , . , , „ , , 171m 

r rm.' m ■ 1 j i. ± , rij. With an estimated density ~ 0.4 pj, WASP-4b ap- 

vest so far. 1 his efficiency is not only due to the constant , , J ; J ' 1 

. . , . ,,, , i j r ii , peared to belong to the subgroup of the planets with a 

optimization of the observational and follow-up strategy, ^ ° , , , . ,f r . ,• , 

, , . ii, i . /r , i onn «N , , radius larger than predicted by basic models of irradiated 

reduction and data analysis (Cameron et al. 2007), but _° ^ V, , 

i , , , c ' , r, , i £ planets (Burrows et al. 2007a, fcortney et al. 2007) unlike 

also to the recent starting of the Southern counterpart of f TTAr „ _\ . , _ . ,„. ' 

the SuperWASP-North facility. Located at the Sutherland ^ASP-Sb (p 7 J, Ji g . Y , 

than WASP-5b (1.2 vs 1.6 Mj), while both planets have a 

similar irradiation (1.89 vs 1.92 10 9 erg s _1 cm -2 ), semi- 
Send offprint requests to: michael.gillon@obs.unige.ch . . ___ _ 4TTN , , , 
* t, , , , „ , i .,, , nnnm . maior axis (0.023 vs 0.027 AU) and host star spectral type 
Based on data collected with the FORS2 imager at the , ' r^/nn 

VLT-UT4 telescope and with the UVES spectrograph at the ^ (j7V VS G4V )- 

VLT-UT2 telescope (Paranal Observatory, ESO, Chile) in the Several hypotheses have been proposed to explain the 
programme 280.C-5003. radius anomaly shown by some highly irradiated planets 
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such as WASP-4b (see Guillot 2008), most importantly 
tides (Bodcnheimer et al. 2001; Jackson et al. 2008), tides 
with atmospheric circulation (Guillot & Showman 2002) 
and enhanced opacities (Guillot et al. 2006, Burrows et 
al 2007a). Receiving similar irradiation from their host 
stars while having significantly different radii, WASP-4b 
and WASP-5b represent a good test for theory and an 
interesting opportunity of progress on our understanding 
of the radius heterogeneity observed among the highly ir- 
radiated planets. It is thus desirable to obtain for these 
two planets the highest precision possible on the system 
parameters. This motivated us to use the VLT to obtain 
(1) a high cadence high precision transit light curve with 
the FORS2 camera, and (2) a high resolution spectrum of 
the host stars with the UVES spectrograph. We present 
respectively in Sec. 2 and 3 these new VLT photometric 
and spectroscopic observations and their reduction. We 
analyze these new data in combination with former transit 
photometry and RV measurements in Sec. 4. The results 
of our analysis are discussed in Sec. 5. 

2. VLT/FORS2 Transit Photometry 

2.1. WASP-4 

The photometry for WASP-4 was obtained on October 23, 
2007. Altogether 339 exposures were acquired with the 
FORS2 camera on the VLT/UT4 telescope from OlhOl 
to 06h50 UT. To have enough reference flux to properly 
correct the photometry from atmospheric effects, the stan- 
dard resolution mode was used, resulting in a 6.8' x 6.8' 
field of view. To obtain a good time sampling of the light 
curve, a 2x2 binning of the pixels was performed, result- 
ing in a pixel scale of 0.25" /pixel. The exposure time was 
tuned by the ESO staff astronomer to 20s while the mean 
read-out plus overhead time was 34s. We chose to observe 
in the z-GUNN+78 filter (X eff = 910 nm, FWHM = 130.5 
nm) to minimize the impact of the stellar limb-darkening 
uncertainty on the deduced system parameters. A very 
large defocus was used to obtain a good duty cycle and 
to minimize the influence of flat-fielding errors: the mean 
characteristic profile width was 50 pixels = 12.5" . Despite 
this large defocus, there was no PSF overlap for the tar- 
get and the reference stars. The guiding system was turned 
on to make the stellar fluxes registered on nearly the same 
pixels during the run (ccntroid jitter ~ 3.5 pixels for the 
whole run) . There was an interruption of 25 minutes in the 
sequence due to a technical problem with the secondary 
mirror setting, fortunately before the transit. The quality 
of the night was photometric. The moon illumination was 
87 %. It was at 44° at closest from the target. The airmass 
decreased from 1.08 to 1.05 then increased to 1.95 during 
the run (Fig. 1). The defocus was tuned several times to 
adapt it to atmospheric transparency variations due to the 
increase of airmass. 

After a standard prc-rcduction, the stellar fluxes were 
extracted for all the images with the IRAF DA0PH0T aper- 
ture photometry software (Stetson, 1987). As the defocus 



was not the same for the whole run, the reduction param- 
eters were adapted to the characteristic profile width of 
each image. 

The transit is already very clear in the absolute flux 
curve of WASP-4. But as shown in Fig. 2, the absolute 
photometry of WASP-4 and of several other stars in the 
field suffer from an unexpected effect: while the shape of 
the largest part of these curves shows a nice airmass-flux 
correlation indicating that the night transparency condi- 
tions were very good, the first part seems to be affected 
by a large systematic dependent on the position on the 
chip. High-accuracy transit photometry has already been 
obtained with the FORS cameras (e.g. Gillon et al. 2007a, 
Pont et al. 2007), and this systematic was not detected in 
these former data. The main difference between these for- 
mer observations and ours is the large defocus that we 
used. The explanation that we and the ESO staff favor 
is linked to the strange shape taken by the primary mir- 
ror Ml with respect to the secondary M2 in case of out- 
of-focus observations, with variations induced by the dif- 
ferent amount of tangential component of the gravity as 
the dominant effect. Indeed, the active optics system of 
the telescope is supposed to compute for each exposure 
an optimal shape for the Ml so to correct for tangential 
gravity pull, and actuators perform micrometrical shifts 
of the Ml to obtain the computed shape, but here the ac- 
tive optics system was turned off at the beginning of the 
run to obtain the required huge defocus. Thus the needed 
correction for the tangential gravity pull was not applied. 
The resulting spatial difference in illumination of the chip 
could then have been rather large at the beginning of the 
run, when the telescope was close to the meridian. As it 
moved away from meridian, the tangential correction be- 
came less important, and so did the effect. 

Differential photometry was performed using the flux 
of several bright stable stars in the field, but the obtained 
curve is still plagued with a large systematic in its first 
part (see Fig. 3). Fortunately, the transit occured in the 
second part of the run for which the effect seems to be 
absent. We thus decided to use only the data after BJD 
= 2454396.625, for which the photometry seems to be re- 
liable and accurate. 

After a careful selection of the reference stars and the 
reduction parameters, we subtracted a linear fit for magni- 
tude vs airmass to correct the photometry for differential 
reddening using the out-of-transit (OOT) data. The cor- 
responding fluxes were then normalized using the OOT 
part of the photometry. The resulting transit light curve 
is shown in Fig. 4, with the best-fit transit model (see Sec. 
4) superimposed. The rms of the first OOT part is 420 
ppm. This value is very close to the theoretical error per 
point obtained from the photon noise of the target and the 
reference stars, the sky background, read-out and scintilla- 
tion noises (Gilliland et al. 1993): 400 ppm. For the second 
OOT part, the measured rms is 740 ppm while the median 
theoretical error is 510 ppm. This largest discrepancy be- 
tween both values comes probably from the amplification 
of the effect of any transparency inhomogeneity across the 
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Fig. 1. Evolution of the airmass during the WASP-4 
VLT/FORS2 run. 
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Fig. 2. Top left : normalized absolute VLT/FORS2 pho- 
tometry for WASP-4. Top right and bottom: normalized 
absolute photometry for other stars of the field. 



field at high airmass. Indeed, the airmass ranges from 1.45 
to 1.95 in the second OOT part. 

The rms and the time sampling are not the only pa- 
rameters needed to evaluate the quality of a photometric 
time series, the level of low-frequency noise (red noise) has 
also to be taken into account, especially for high SNRs 
(Pont et al. 2006). We estimated the level of red noise ay 
in our photometry using the equation (Gillon et al. 2006): 



I N-l 



1/2 



(1) 



where a is the rms in the original OOT data and o~n is 
the standard deviation after binning the OOT data into 
groups of N points. We used N = 20, corresponding to a 
bin duration similar to the ingress/egress timescale. The 
obtained value for a r is quite small: 110 ppm. 
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Fig. 3. Differential photometry for WASP-4 before rejec- 
tion of the part of the curve damaged by the systematic 
and differential reddening correction. 



WASP4 VLT/F0RS2 z-band 




o 

-0.002 



0.001 

o 

-0.001 



1 , . , . 1 , . , . 1 


4396.65 


BJD 


4396.7 
- 2454000 (days) 


4396.75 






L ' ' | ' 






'o 1 1 






'- . , i 




i 


o o 

1 





o 


4396.65 


BJD 


4396.7 
- 2454000 (days) 


4396.75 








i j 


{ I 1 








= , , i 




i 


i 


i 




4396.65 


BJD 


4396.7 
- 2454000 (days) 


4396.75 







Fig. 4. Top: VLT/FORS2 z-band transit photometry for 
WASP-4. The best-fit transit curve is superimposed in red. 
Middle: residuals of the fit (rms = 570 ppm ). Bottom: 
residuals of the fit after binning per 20 points (rms = 200 
ppm) . 

2.2. WASP-5 

The photometry for WASP-5 was obtained on November 
15, 2007. 337 exposures were acquired from 00h24 to 05h45 
UT, again with the FORS2 camera. The same observa- 
tional strategy as for WASP-4 was used. The quality of the 
night was photometric. The moon was in its first quarter 
(illumination = 24 %), and its closest distance to the tar- 
get was 48° at the end of the run. The airmass decreases 
from 1.06 to 1.04 then increased to 1.97 (Fig. 5). 

The same reduction procedure as for WASP-4 was 
used. Here again, the photometry of the first part of the 
run shows the presence of a systematic at the percent 
level (Fig. 6). We notice that the amplitude of the ef- 
fect depends here too on the position on the chip and ap- 



4 



M. Gillon et al.: Improved parameters for WASP-4b and WASP-5b 




419.S5 



419.6 419.65 419.7 

BJD - 2454000 (days) 



Fig. 5. Evolution of the airmass during the WASP-5 
VLT/FORS2 run. 
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Fig. 6. Top left : normalized absolute VLT/FORS2 pho- 
tometry for WASP-5. Top right and bottom: normalized 
absolute photometry for other stars in the field. 



pears to affect only the low-airmass photometry, it is thus 
probably the same effect than for WASP-4 (see above). 
Unfortunately, the transit happened at the beginning of 
the run and its photometry is too damaged to be useful. 
We thus decided to use only the UVES spectroscopy (see 
Sec. 3) and former data to characterize this system. 



3. VLT/UVES Spectroscopy 

High-resolution spectra of WASP-4 and WASP-5 were ob- 
tained using the UVES spectrograph on the VLT/UT2 
telescope on November 27, 2007. The red arm was used 
with standard setting 580, giving spectral coverage from 
4780 to 6808 A, except from 5758 to 5833 A due to the gap 
between to two CCD detectors. A 1-arcsec slit was used 
giving a spectral resolution of 40 000. Exposure times of 
600 s were used yielding mean S/N of 94:1 and 102:1 for 



WASP-4 and WASP-5, respectively. The standard pipeline 
reduction products were used in the analysis. 

The analysis was performed using the uclsyn spectral 
synthesis package (Smith 1992; Smalley et al. 2001) and 
ATLAS9 models without convective overshooting (Castelli 
et al. 1997). The H a line was used to determine the effec- 
tive temperature (T e g), while the Na I D and Mg I lines 
were used as surface gravity (log g) diagnostics. The pa- 
rameters obtained from the analysis are listed in Table 1. 

The equivalent widths of several clean and unblended 
lines were measured. Atomic line data was mainly taken 
from the Kurucz & Bell (1995) compilation, but with up- 
dated van der Waals broadening coefficients for lines in 
Barklem et al. (2000) and log gf values from Gonzalez 
& Laws (2000), Gonzalez et al. (2001) or Santos et al. 
(2004). A value for microturbulence (£ t ) was determined 
from Fe I using Magain's (1984) method. The ionization 
balance between Fe I and Fe n and the null-dependence of 
abundance on excitation potential were used as additional 
T e g and logg diagnostics (Smalley 2005). 

In addition to the spectral analysis, we have also used 
published broad-band photometry to estimate the total 
observed bolomctric flux (/©). For WASP-4, TYCHO- 
2 (H0g et al. 2000), USNO-B1.0 R-mag. (Monet et al. 
2003) and 2MASS (Cutri et al. 2003) were used, while 
for WASP-5 GALEX nuv flux (Morrissey et al. 2007), 
NOMAD (Zacharias et al. 2004), DENIS (Fouque et al. 
2000) and 2MASS were used, but no TYCHO photome- 
try was available. The photometry was converted to fluxes 
and the best-fitting Kurucz (1993) model flux distribution 
found, which was integrated to determine /®. The Infrared 
Flux Method (Blackwell & Shallis 1977) was then used 
with 2MASS magnitudes to determine T c g and stellar an- 
gular diameter (9). The results are given in Table 1 and 
are consistent with those from the spectroscopic analysis. 

We have determined the elemental abundances of sev- 
eral elements (listed in Table 1) from their measured 
equivalent widths. The quoted error estimates include 
those given by the uncertainties in T e ff, log g and £ t , as 
well as the scatter due to measurement and atomic data 
uncertainties. In our spectra the Li I 6708A line is not de- 
tected. Thus we can only give upper-limits on the Lithium 
abundances. 



4. Data Analysis 

4.1. Determination of the system parameters 
4.1.1. WASP-4 

We derived stellar and planetary parameters for the sys- 
tem by fitting simultaneously our new VLT z-band transit 
light curve with the data presented in W08, i.e. (1) the 
WASP R-band photometry, (2) the FTS i-band transit 
light curve, (3) an EulerCAM R-band transit light curve 
and (4) 14 CORALIE radial velocity measurements. 
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The data were used as input into the Monte Carlo 
Markov Chain (MCMC) code described in Cameron et al. 
(2007) which was designed specifically to solve the multi- 
variate problem of transiting star-planet systems. Via the 
MCMC approach, the fitting code repeatedly adopts trial 
parameters until it converges on a set of values which 
produces the best model velocity curve and model light 
curves. In short, nine parameters were used to describe 
the light curves and radial velocity curve of the host star 
including the orbital period P, the time of minimum light 
T , the transit depth 5, the total transit duration tr, the 
impact parameter b, the stellar mass M*, the stellar ve- 
locity amplitude K\, the systemic radial velocity 7, the 
orbital eccentricity e, and the longitude of periastron w. 
These nine fitted parameters determine the physical prop- 
erties of the star-planet system, including the masses and 
radii of the star and planet and the orbital inclination 
and separation. The goodness-of-fit statistic used to assess 
the best parameters is the sum of the x 2 for all the data 
curves with respect to the models. Model light curves were 
derived according to the formalism outlined in Mandel 
& Agol (2002), adopting the small-planet approximation 
and using the non-linear limb darkening coefficients from 
Claret (2000, 2004) for the appropriate photometric fil- 
ters. The code also has the option to apply a Bayesian 
main sequence prior on the stellar mass and radius which 
acts to keep the star on the main sequence. However, due 
to the exceptional quality of the follow-up photometry, we 
did not apply this constraint in the analysis of WASP-4. 

We ran the MCMC code in an iterative fashion to de- 
rive the best overall solution for the properties of the star- 
planet system. We also combined the physical properties 
derived from the MCMC code with the stellar parameters 
from the spectral synthesis to determine the evolutionary 
status of the host star and to confirm that the stellar prop- 
erties, including mass (which is not directly measured), 
were consistent with each other and with theoretical stel- 
lar evolution models. 

In the initial run of the MCMC code, we adopted ini- 
tial guesses for the light curve parameters from the results 
of the box-least squares analysis of the WASP data. We 
also assumed a starting value for the eccentricity of 0.02, a 
systemic RV equal to the mean of the velocity data, and a 
velocity amplitude derived by fitting a sinusoidal velocity 
variation to the observed RVs by minimizing \ 2 ■ The ini- 
tial guess for the stellar mass was derived by interpolating 
the zero age main sequence, solar metallicity stellar evolu- 
tion isochrone of Girardi et al. (2000) at the temperature 
of the host star. In addition, the input light curves had the- 
oretically derived uncertainty values for each photometric 
measurement. The uncertainties were computed by con- 
sidering shot noise, scintillation, read-out and background 
noises, but did not include correlated noise. For the RV 
measurements, the MCMC code adds quadratically an ad- 
ditional RV jitter to the theoretical uncertainties so that 
the reduced \ 2 °f the RV data compared to the model 
curve is approximately equal to 1 . Finally, we allowed the 
eccentricity to float freely. 



The results of this initial run were used to inform the 
second and final run of the code in the following ways. 
First, the initial eccentricity result was well within la of 
zero. Thus, for final run, we fixed the eccentricity to zero, 
assuming a circular orbit for the planet. Next, the best 
fitting model light curves from the initial run were used to 
make a correlated noise measurement for each photometric 
time-series. We employed the method described in Sec. 2.1 
(Eq. 1), but using the entire residual time-series, rather 
than just the OOT parts. Table 2 presents the red noise 
values obtained. For each photometric measurement, we 
added the appropriate correlated noise value in quadrature 
to the theoretical uncertainty. The resulting light curves 
were input into the final run of the MCMC code. 

Finally, the first MCMC run produced a measurement 
of the mean stellar density (p s = M s /E? s ) which we used 
to refine the initial guess of the stellar mass. We converted 

1 /3 

the derived stellar density to R s /M s in solar units, and 
compared this property and the stellar temperature in a 

modified Hertzsprung-Russel (HR) to the Girardi solar 

1/3 

metallicity models. The quantity, R s /M s ' , depends only 
on the observed transit properties (duration, depth, im- 
pact parameter, and orbital period) and is independent of 
the measured temperature. We generated the same prop- 
erty from the mass and log g values in the models, and 
then interpolated the models in the i?/M 1 / 3 -Teff plane to 
determine a mass and age for WASP-4. We interpolated 
linearly along two consecutive mass tracks to generate an 
equal number of age points between the zero-age main se- 
quence and the evolutionary state defined as Te-M which 
is the stage where the star reaches core Hydrogen exhaus- 
tion. We then interpolated between the mass tracks along 
equivalent evolutionary points to find the mass and age 
from the models that best match the stellar properties 
derived from the MCMC code (density) and the spectral 
synthesis (temperature). In this way, we obtained a new 
intital guess for the stellar mass of, M* = 0.93 M Q . 

We then implemented a second and final run of the 
MCMC code, (1) fixing e = 0, (2) including correlated 
noise in the photometric uncertainties, (3) applying the 
results from the first MCMC run as the initial guesses 
for the fitted parameters, and (4) adopting a prior on the 
stellar mass of M* = 0.93 M@. The results of this run are 
given in Table 3. The obtained value for the RV jitter is 7 
m s _1 , identical to the value deduced in W08. 

Lastly, we plotted the final stellar parameters on the 
modified HR diagram and interpolated the model tracks 
as described above to determine the age of the system. We 
derive an age of 5.2+32 Gyr. We also confirmed that the 
final mass derived in the MCMC code was consistent with 
the observed temperature and the stellar evolution mod- 
els (Figure 7). The MCMC result gives a mass of 0.85^;^ 
Mq which is within the la uncertainty on the mass de- 
termined from the theoretical tracks (0.93 ± O.O5M ). 

We also made an independent analysis aiming to test 
if the quality of our VLT z-band transit photometry is 
good enough to constraint reliably the limb-darkening of 
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units versus effective temperature. The properties of the 
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are labeled and the isochrones are 100 Myr (solid), 1 Gyr 
(dotted), 5 Gyr (dot-dashed), 10 Gyr (dotted). According 
to the models, the host star has an age of 5.2±|;| Gyr. 
The y-axis error bars are smaller than the data point. 



the star. In this analysis, we used as data our VLT tran- 
sit light curve with only the WASP photometry and the 
CORALIE radial velocities. We assumed a quadratic limb- 
darkening law for the photometric models to minimize 
the number of free parameters and allowed the two limb- 
darkening coefficients u\ and U2 to float for the VLT pho- 
tometry only, using as jump parameters not these coeffi- 
cients themselves but the combinations 2 x u\ + u 2 and 
u\ — 2 x ii2 to ensure that the obtained uncertainties arc un- 
correlated (Holman et al. 2006). We obtain values for the 
transit depth and duration that are slightly different from 
the values shown in Table 3, but all the deduced physi- 
cal parameters are in good agreement with the previous 
values. The obtained values for the limb-darkening coef- 
ficients are m = 0.299^^26 and u 2 = 0.2481^23- These 
values are physically plausible in the sense that they pro- 
duce a monotonically decreasing intensity from the center 
of the star to the limb, but they are consistent with the 
values interpolated from Claret's tables (2000, 2004): u\ 
= 0.266 and u 2 = 0.302. This disagreement could be due 
to the fact that our VLT photometry has a small but sig- 
nificant level of correlated noise able to modify slightly the 
actual shape of the transit. The amplitude of this red noise 
is in fact similar to the difference between the models fit- 
ted with and without free limb-darkening coefficients. We 
thus prefer to consider the values for the system param- 
eters obtained with a fixed non-linear limb-darkening law 
as our final ones. 



4.1.2. WASP-5 

An analysis similar to the one described above was per- 
formed for WASP-5 using data presented in A08: (1) the 
WASP photometry, (2) a FTS i-band transit light curve, 
(3) an EulerCAM R-band transit light curve and (4) 11 
CORALIE radial velocities. As in A08, no additional RV 
jitter was needed to obtain a final reduced \ 2 close to 1 
for the RV model. No Bayesian main sequence prior was 
used as for WASP-4. The orbital eccentricity obtained in 
the first MCMC run was non-zero at the 2a level, so it 
was kept as free parameter in the second run, leading to 
a final value of 0.04 ± 0.02 (see Table 3). Nevertheless, e 
and oj are not orthogonal parameters because e cannot be 
< 0, and this asymmetry could lead to an overestimation 
of e for low eccentricity orbit. To check that our marginal 
detection of a non-null eccentricity is reliable, we made a 
new analysis using as jump parameters e.cosw and e.sinu>, 
which are orthogonal, and calculated e and u> afterwards. 
We obtained e = 0.049^;^ and uj = 0.73±°;^, in very 
good agreement with the values presented in Table 3. We 
thus claim the marginal detection of a non-null eccentric- 
ity for WASP-5b. 

WASP-5 is slightly super solar metallicity with 
[Fe/H]= +0.09 ± 0.09 (see Table 1), so mass tracks of the 
Girardi models in metallicity at +0.09 were first linearly 
interpolated between the zero and +0.20 metallicity mass 
tracks before the obtained stellar parameters were plotted 
on the modified HR diagram and model tracks were inter- 
polated (Fig. 8). From this final step, we derive an age of 
5. 4+43 Gyr and a stellar mass of l.OO±8;8£M , within the 
error bar of the MCMC value of 0.96±g"J| M Q (see Table 
3). 

4.2. Transit timings 

Fixing all the system parameters except epoch to the 
ones deduced from the above analysis, we fitted with our 
MCMC code a transit profile to each transit to obtain in- 
dividual timings. For the WASP data, the transits of the 
season 2006 and 2007 were folded together to obtain a tim- 
ing per season. Table 2 gives the obtained timings and er- 
ror bars. From these results, the VLT z-band photometry 
for WASP-4 appears to show a formal error of <~ 5s on its 
deduced timing, comparable to the best precision obtained 
from space (see e.g. Knutson et al. 2007). Such a precision 
is doubtful for our ground-based data, because we know 
that despite their high quality they have a low but still 
significant level of covariant noise able to bring a system- 
atic error on the deduced timing. To assess the influence of 
this covariant noise on our timing precision, we analyzed 
each transit with the method described in Gillon et al. 
(2007b), for which the estimation of the errors is based on 
the 'prayer bead' procedure: after having determined the 
best-fitting eclipse model, a large number of fits are per- 
formed and for each of them the residuals of the initial fit 
are shifted sequentially about a random number and then 
added to the eclipse solution. This procedure allows to 
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Fig. 8. Similar modified HR diagram than Fig. 7 but here ^ o ' 
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take into account the actual covariant noise of the data. At 
the end, the error bars of the fitted parameters are deter- 
mined from the distribution of their derived values. Table 
2 presents also the values and error bars obtained with this 
method. We notice that for the data showing the largest 
level of covariant noise, the WASP data, the error bars on 
the timing is ^4 times larger. For the Euler transits having 
a rather low level of covariant noise, the difference is much 
smaller. Interestingly, the error bars on the WASP-4 VLT 
transit timing is 3 - 4 times larger than the one obtained 
with the MCMC code and is only slightly better than the 
one of the Euler transits. This shows clearly that covari- 
ant noise has an important impact on the transit timing 
precision. While Table 2 presents the timings and errors 
obtained with both methods, we outline that our final re- 
sults are the ones obtained with the 'prayer bead' method: 
only these values should be used in future analysis. 

Figure 9 shows for both planets the residuals from the 
subtraction to these timings of the calculated transit tim- 
ings based on ephemeris presented in Table 3. The transits 
of WASP-4 show no clear sign of period variability. A lin- 
ear fit to the transit timings as a function of the transit 
epoch results in a period of 1.3382319 ± 0.0000068 days, 
in excellent agreement with the period obtained from our 
combined MCMC analysis. The reduced x 2 01 the fit is 
1.35, indicating that the epoch - timing relation is well 
modeled by a line. For WASP-5, the period deduced from 
a similar fit is 1.628430 ± 0.000013 days, in good agree- 
ment with the value obtained from the combined analysis, 
but here the reduced x 2 is 5.7. At this stage, we cannot 
assign any firm significance to a possible period variabil- 
ity, because a shift of one of the two most precise timings 
(FTS or Euler) due to an unknown systematic effect could 
lead to a much bettered reduced x 2 , but obtaining more 
precise transit timings for WASP-5 is desirable. 



Fig. 9. Observed minus calculated (O-C) transit timings 
for the light curves included in this study. Table 2 lists 
the actual transit timings. The calculated timings were 
obtained from the ephemeris due to our combined MCMC 
analysis (see Table 3). 

5. Discussion 

5.1. Ground-based photometry 

Until recently, it was considered by many that ground- 
based photometry could not reach the high cadence sub- 
mmag regime because of the presence of the atmosphere. 
Indeed, high frequency atmospheric noises (mainly scintil- 
lation) limit the precision that high SNR photometry can 
reach within small time bins. If one is willing to compro- 
mise on the sampling of their photometric time-series, bin- 
ning the data (or using longer exposures) allows for getting 
better errors, but the obtained precision will be finally lim- 
ited by low frequency noises. To observe several times the 
same planetary eclipse and to fold the photometry with 
the orbital period is thus generally considered as the only 
option to get very well sampled and precise eclipse light 
curves from the ground. Nevertheless, we show here that 
reaching the sub-mmag sub-min regime for one eclipse is 
possible with a large aperture ground-based instrument. 
Photon noise and scintillation are not a concern even for 
small time bins when a rather bright transiting system like 
WASP-4 is monitored with a large aperture telescope like 
the VLT. The high standard quality of such a telescope 
and the excellent atmospheric conditions at Paranal lead 
furthermore to a very low level of covariant noise in the 
differential photometry. 

Unfortunately, we report the presence of an instru- 
mental effect on the VLT that damaged a part of our 
photometry and that could be a major problem for simi- 
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lar programs. It was unfortunately the case for WASP-5: 
the photometry is too damaged by the effect to be useful. 
As mentioned in Sec, 2, the effect is probably due to the 
fact that we had to turn off the active optics system of the 
VLT to obtain the required very large defocus. We empha- 
size that imaging in this manner is not at all a standard 
observational mode on the VLT. We suggest for similar 
programs a milder level of defocus combined with the use 
of the active optics system. 

Recently, exquisite ground-based transit photometry 
for the hot Neptune GJ 436b was obtained by Alonso et al. 
(2008) using a different approach. They observed a tran- 
sit in the H-band with the TCS telescope and its CAIN-II 
near-IR detector. As the red dwarf GJ 436 is very bright in 
the H-band (H = 6.3), the background variability is not a 
concern, so no dithering pattern was used and the images 
were severely defocused, i.e. the strategy was very similar 
to the one that we choose for our VLT observations. The 
most surprising point is that no differential photometry 
was used by Alonso et al. to reach such a high photo- 
metric quality. These authors explain this by the much 
smoother behavior of the transparency variations in the 
H-band compared to the visible. It is very desirable to 
confirm this point by obtaining more high-quality eclipse 
light curves in the near-IR. Unfortunately, this method is 
limited to stars that are very bright in the near-IR, and 
only a few are known to harbor a transiting planet (e.g. 
HD 189733, HD 209458). 

5.2. The hot Jupiters WASP-4b and WASP-5b 

The high-quality of the WASP-4 VLT z-band transit pho- 
tometry allows a significant improvement of the precision 
on the impact parameter (0.063±^ vs 0.13±g;^ in W08) 
that leads to a better precision on the orbital inclination 
(89.351^9 vs 88.59±^ degrees in W08), but our anal- 
ysis fails to give a more precise value for the planetary 
mass (1.21+g;£f Mj vs 1.215±g;g| Mj in W08) and even 
for the planetary radius (1.304±g;g|^ Rj vs 1.416±g;g2i 
Rj in W08) despite a significantly more precise determi- 
nation of the planet/star area ratio (0.02357i. 00O0O8 vs 
0.01192±g;ggg|f Rj in W08). This is due to the fact that 
our new spectroscopy docs not improve significantly our 
knowledge of the host star because it is now limited by 
the accuracy /validity of the stellar atmospheric and evo- 
lution models, and this uncertainty on the stellar param- 
eters propagates to our final accuracy on the planet pa- 
rameters. This fact is also illustrated with our results for 
WASP-5b: our new values for the planetary mass and ra- 
dius agree well with the ones quoted in A08 (1.58+ ; 10 Mj 
and 1.087±g;gff Rj vs 1.58±g;£| M J and l-090±g;g§| Rj 
in A08) but are not more accurate. This 'stellar' accuracy 
limit relates to all transiting systems, except the putative 
systems (1) that could be studied by asteroseismology (see 
e.g. Christensen-Dalsgaard et al. 2007), allowing a more 
precise estimation of their age and radius, or (2) the sys- 
tems that would be bright and nearby enough to allow a 



very precise direct determination of the stellar radius by 
long baseline interferometry (see e.g. Baines et al. 2008). 

With a radius measurement R p = 1.301q'q4 we 
confirm here that WASP-4b is larger than predicted by 
basic models of irradiated planets (Burrows et al. 2007a, 
Fortney et al. 2007). For instance, the theoretical value 
presented by Fortney et al. (2007) for an irradiated 1.46 
Mj core-less planet of 4.5 Gyr orbiting a 0.02 AU from a 
sun- like star is only 1.17 Rj. WASP-4b is thus another 
case demonstrating that something is missing in basic 
models. 

Fortney et al. (2008) proposed the theoretical divi- 
sion of hot Jupiters into two classes based on their level 
of irradiation. By analogy with M-dwarfs, the 'pM' class 
would be composed of the planets warmer than required 
for condensation of Ti and V-bearing compounds, and 
these planets should show a stratospheric temperature in- 
version due to the absorption of most of the large inci- 
dent flux by the high-opacity gaseous TiO and VO. The 
cooler planets would compose the 'pL' class. Burrows et 
al. (2008) proposed a similar bifurcation into two groups, 
but without firmly identifying TiO and VO as the high 
opacity gaseous compounds causing the thermal inversion. 
Interestingly, recent Spitzer secondary eclipse measure- 
ments detected such temperature inversions for the highly 
irradiated planets HD 209458b and XO-lb (Burrows et al. 
2007b, Knutson et al. 2008, Machalek et a. 2008). The up- 
dated estimation for the incident flux received by WASP- 
4b is now 1.7 10 9 erg s^ 1 cm~ 2 , classing it clearly in the 
theoretical pM planetary class proposed by Fortney et al. 
Figure 10 shows the location of WASP-4b in a mass-radius 
diagram. The location of the other transiting planets re- 
ceiving an incident flux larger than 10 9 erg s _1 cm~ 2 (i.e. 
belonging clearly to Fortney's pM class) or smaller than 5 
x 10 8 erg s _1 cm~ 2 (i.e. belonging clearly to Fortney's pL 
class) is also shown for comparison. We notice that the pL 
planets seem to have a smaller radius than their pM coun- 
terparts having a similar mass. WASP-4b follows this ten- 
dency. This is in favor of the proposition done by Burrows 
et al. (2007a) that the enhanced opacity of the highly irra- 
diated giant planets could alter their cooling history and 
thus be part of the solution to their 'anomalously' large 
radius. However, WASP-5b is even more irradiated then 
WASP-4b (2.0 x 10 9 erg s~ x cm~ 2 ) but according to our 
result it is not extremely bloated (R p = 1.09 ± 0.07 Rj). 
As shown in Fig. 10, WASP-5b seems to be smaller than 
the other planets of similar mass that fall clearly in the 
pM class. This favors the fact that the size of highly irra- 
diated gazeous planets is not only dependent on their level 
of stellar irradiation but also of other factors. The plane- 
tary core mass is probably one of them. In this context, 
it is interesting to notice that WASP-5 is more metal-rich 
than WASP-4. This goes in the right direction towards 
the existence of a correlation between the metallicity of 
the parent star and the core mass of the planet (Guillot 
et al. 2006, Burrows et al. 2007a, Guillot 2008). 

More transit photometry for WASP-5 is needed to con- 
firm the possible period variability shown by our transit 
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Fig. 10. Location of WASP-4b (W4) and WASP-5b (W5) 
in a mass-radius diagram. The location of the other tran- 
siting planets receiving an incident flux larger than 10 9 
erg s" 1 cm" 2 (open circles) or smaller than 5 x 10 8 erg 
s _1 cm" 2 (closed circles) is also shown. 

timing measurements. Such an apparent period variabil- 
ity could be due to the presence of another close-in lighter 
planet in the system (see e.g. Holman & Murray 2005). 
A precise timing of the secondary eclipse of this planet is 
also desirable, because it could confirm the non-null ec- 
centricity that we marginally deduce from our analysis. 
Such a non-null eccentricity could also be the sign of the 
presence of another planet around WASP- 5. 
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WASP-4 


WASP-5 


Toff- 


5500 ± 100 K 


5700 ± 100 K 


logS 


4.5 ± 0.2 


4.5 ± 0.2 


& 


1.1 ± 0.2 km s" 1 


1.2 ± 0.2 km s" 1 


v sin i 


2.0 ± 1.0 km s" 1 


3.5 ± 1.0 km s" 1 


[Fe/H] 


-0.03 ± 0.09 


+0.09 ± 0.09 


[Si/H] 


+0.08 ± 0.05 


+0.21 ± 0.06 


[Ca/H] 


+0.04 ± 0.14 


+0.07 ± 0.14 


[Sc/H] 


+0.07 ± 0.13 


+0.21 ± 0.10 


[Ti/H] 


+0.11 ± 0.12 


+0.14 ± 0.08 


[V/H] 


+0.09 ± 0.07 


+0.16 ± 0.07 


[Cr/H] 


+0.06 ± 0.08 


+0.07 ± 0.10 


[Co/H] 


+0.07 ± 0.07 


+0.20 ± 0.07 


[Ni/H] 


+0.02 ± 0.10 


+0.15 ± 0.05 


log N(Li) 


< 0.8 


< 0.5 


Teff(IRFM) 


5470 ± 130 K 


5740 ± 130 K 


0(IRFM) 


0.031 ± 0.002 mas 


0.032 ± 0.002 mas 



Table 1. Spectroscopic parameters obtained in this work for WASP-4 and WASP-5. 



Light curve 



oy [ppm] 



Mid-transit timing (MCMC) [BJD] 

"~ " H). 00074 



Mid-transit timing ('prayer bead' 
_ .,+0.0025 



[BJD] 



WASP-4 WASP06 R-band 3150 

WASP-4 WASP07 R-band 4020 

WASP-4 Euler R-band 

WASP-4 FTS i-band 510 

WASP-4 VLT z-band 190 



2453963.10863I 0(VI| . 

2454364.57722t°;°g^ 
2454368.59244t°'°°^2 

2454371.26812±°;°°°g 
2454396.695410 ± 0.000051 



2453963.1086 
2454364.5757 
2454368.59266 
2454371.26738 
2454396.69548 



0.0023 
+0.0021 
0.0033 
+0.00025 
0.00027 
0.00097 
0.00087 
+0.00015 
0.00026 



WASP-5 WASP06 R-band 3260 

WASP-5 WASP07 R-band 1490 

WASP-5 Euler R-band 360 

WASP-5 FTS i-band 790 



2453945.71962+^;™™^ 

2454364.2283tg ; ^3 
2454383.76684t°'°°°25 



2454387.02221 



+0.00034 



2453945.7187Z 
2454364.2285 

2454383.76738; 

2454387.02197 



+0.0057 
-0.0069 
-0.00031 
-0.00032 
+0.00071 
-0.00050 



Table 2. Typical red noise values and mid-transit times evaluated for the photometric time-series used in this analysis. 



Parameter 



Symbol 



WASP-4 



WASP-5 



2454383.313070t£S£ 

1 qqooqo/i +0.0000017 
l.OOOiO/;i_ .0000029 

fl ri9o.c;7+ ' (,()()1() 
U.UZOd I _().00008 

0.088311°:°™? 
0.063±g;g£ 



2454373.99598loSj;jiq 

1 R9S/L97Q+0. 0000022 
1 ' uimil a -0. 0000049 

o.oii8ol° : °r29 

0.09871°:°^ 



Units 



Transit epoch (HJD) 
Orbital period 
Planet/star area ratio 
Transit duration 
Impact parameter 



P 

(R P /Rs) 2 
t T 

b 



0.31 



+0.06 



days 
days 

days 



Stellar reflex velocity 
Centre-of-mass velocity 
Orbital semi-major axis 
Orbital inclination 
Orbital eccentricity 
Longitude of periastron 



7 
a 
/ 

e 



2476+ ' 0139 

"•^'"-O.OOHS 
c 7 V qo 7 +0.0020 
01 -lOOl _0.0013 

fl n99c;c;+ (l ' ()(,(l!,r> 
U.UZZdO_ .00065 

89.35lg:S 
0.0 (fixed) 



9707+0.0092 

U.z f »( -0.0063 


km s 1 


20.0087l°;™ 3 25 


km s _1 


n267 +0 '°° 12 

U.UZO( _o.0008 


AU 


86.9+07 


degrees 


0.0381°:'°?? 




o.6olg;S 


rad 



Stellar mass 
Stellar radius 
Stellar surface gravity 
Stellar density 



R* 
log g* 

P* 



0.85 



+0.11 
-0.07 
+0.036 



0. 873_o.o27 
4 487 +o.oi9 

1. tOl _ 015 
1 9S/I+0.013 
J — ao ^ : — 0.019 



0.96 



+ 0.13 



1.029] 



J -0.09 
-, + 0.056 

^-o.oeo 
4 qq^+0-043 

l.O3U_ 040 

0.88 + 0.12 



M 

Rq 

[cgs] 

PQ 



Planet radius 
Planet mass 

Planetary surface gravity 

Planet density 

Planet temperature (A = 0, 



F=l) 



Rp 
M p 
log5p 
Pv 
Teq 



1 qn4+° 054 

l.OUI-o 042 
-I 91 +0.13 
1 - zi -0.08 
q o-i 0+0.025 
O.Z1Z_ on 

S46+ 039 

U.O4O_ 025 

1650 ± 30 



t.US/'-o 071 

1 58+° 13 
i.oo_ io 

o 48 c+0.054 
3.<iOO_ 043 
-. 9 q+0.26 
liZO -0 .16 

1706l|g 



Rj 
Mj 
[cgs] 

pjup 

K 



Table 3. WASP-4 and WASP-5 system parameters and 1-a error limits derived from MCMC analysis. 



